A generic method to constrain the dark matter model parameters from Fermi 

observations of dwarf spheroids 
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Observation of 7-rays from dwarf galaxies is an effective way to search for particle dark matter. 
Using 4-year data of Fermi-LAT observations on a series of Milky Way satellites, we develop a 
general way to search for the signals from dark matter annihilation in such objects. Instead of 
giving prior information about the energy spectrum of dark matter annihilation, we bin the Fermi- 
LAT data into several energy bins and build a likelihood map in the "energy bin - flux" plane. 
The final likelihood of any spectrum can be easily derived through combining the likelihood of all 
the energy bins. It gives consistent result with that directly calculated using the Fermi Scientific 
Tool. This method is very efficient for the study of any specific dark matter models with 7-rays. 
We use the new likelihood map with Fermi-LAT 4 year data to fit the parameter space in three 
representative dark matter models; i) toy dark matter model, ii) effective dark matter operators, 
and iii) supersymmetric neutralino dark matter. 

PACS numbers: 95.35.-|-d,95.85.Pw 



I. INTRODUCTION 

Detecting partiele dark matter (DM) is one of the 
most important task for the modern physics, although 
its existence has been revealed through astronomical ob- 
servations for more than seventy years. The DM par- 
ticles could annihilate themselves to produce standard 
model (SM) particles such as electrons/positrons, pro- 
ton/antiprotons, 7-ray photons and neutrinos, which pro- 
vides us a chance to detect the DM indirectly in the cos- 
mic radiation. The Fermi Large Area Telescope (Fermi- 
LAT) is currently the best detector for 7-ray detection 
in space and could substantially increase the sensitivity 
for DM search. 

The dwarf spheroidal satellites (dSphs) in the Milky 
Way are generally regarded as ideal candidates when 
searching for possible DM annihilation signals due to 
their high DM content and low baryonic contamination. 
The search for 7-ray emission from dSphs with Fermi- 
LAT data was performed in many works With- 
out finding any significant "signal" from these targets, 
stringent upper limits on DM annihilation cross section 
can be derived. It was shown that for < 20 GeV 
the canonical thermally produced DM with cross section 
~ 3 X 10~^^ cm^ s~^ was ruled out with two-year Fermi- 
LAT observations on 10 dSphs ^ (see also [J for similar 
conclusion) . 

There are usually two ways to extract the limits on 
DM models from the 7-ray observations. The most con- 



servative and DM-model-independent way is to require 
the DM- induced signal not to exceed the flux/limit of 
specific source. But in this way, the spectral and spa- 
tial information of DM is missing and the constraint is 
usually weak^ . Another way which could be more realis- 
tic is to incorporate the expected 7-ray spectrum and/or 
the sky map in the analysis. However, it needs to input 
the model of DM when analyzing the data and is not 
convenient for generic purpose. 



In this work we alternatively adopt a new approach to 
extract the limits on the 7-ray emission from the dSphs 
independent with the speculated 7-ray spectrum. The 
Fermi-LAT observational data are binned into a series 
of energy bins, and the likelihood in each energy bin is 
calculated assuming a constant value of the flux from 
DM annihilation in this energy bin. Then the total like- 
lihood of any spectrum can be easily derived with such 
a likelihood map on "energy-flux" plane. The detailed 
description of the method is given in Sec. II. In Sec. Ill 
we apply this method to several DM models. 



^ Note in Q a stronger constraint was given. In that work a lower 
upper limit of 7-ray flux was derived through subtracting the 
photon counts of choosen "background" region from the "signal" 
region. 
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II. LIKELIHOOD MAP 

In this section, we will describe the methodology of 
how to build a particle mo del- independent likelihood 
map. The 7-ray flux from annihilation of DM in a dSph 
is 



HE) 



IK 



X J, 



(1) 



where is the mass of DM particle, (av) is the average 
velocity-weighted annihilation cross section, 



is the 



7-ray yield spectrum of one annihilation of a d\i pair, 
and J = J dldQ.p^{l) is the integral density square of 
DM. 

From Fermi-LAT 2-year result 0, we can sec {av) 
with the bb channel for < 25 GeV is lower than 
the result expected from thermal equilibrium, {av) ^ 
3 X 10^2^ cm^s"^ It could be a hint that Fermi-LAT 
can rule out DM annihilation to bb, at least, at the lower 
mass region, < 25 GeV. However, the above infor- 
mation cannot be true for realistic DM models. In most 
of particle models such as supersymmetry (SUSY), 7- 
rays can be produced via several annihilation channels. 
Moreover, can differ from one point to another in the 
SUSY parameter space. Therefore it will be more useful 
to have a likelihood map which does not depend on a 
pre-determined To do so, we bin the photons into 

small energy bins [Ei,Ei^i], and assume the differential 
energy spectrum in the small energy interval to be 
constant C,. Then the likelihood is built based on two 



new variables, Ei and 



- X d X Jj , where i is 



the index of energy bin and j is the index of individual 
dSph. Wc adopt the same 10 dSphs in this study as in 



To calculate the likelihood from Fermi-LAT data we 
employ 4-year Fermi-LAT data^ recorded from 4 Au- 
gust 2008 to 2 August 2012, with the pass 7 photon 
selection. The LAT Scientific Tools version v9r27pl 
are used for the analysis. The "SOURCE" (evclass=2) 
event class is selected, and the recommended filter 
cut "(DATA_QUAL==1) && (LAT.C0NFIG==1) && 
ABS(ROCK.ANGLE)< 52" is applied. The photon en- 
ergy range is choosen from 200 MeV to 200 GeV, and 
the rcgion-of-intercst (ROI) is adopted to be a 14° x 14° 
box centered around the center of each dSph. The instru- 
ment response function used is "P7S0URCE_V6" . For the 
diffuse background, we use the Galactic diffuse model 
gal_2yearp7v6_v0 . f its and the isotropic diffuse spec- 
trum iso_p7v6source.txt provided by the Fermi Sci- 



Likelihood Map 




FIG. 1: The likelihood map on (Ej, E^ip) plane based on 
4-year Fermi-LAT data on the 10 dSphs. The black solid line 
is one-sided 95% exclusion limit. 



ence Support Center'^. The point sources of the second 
LAT source catalog are also included in the modeling . 
The binned likelihood method is used to calculate the 

likelihood L]"!^'^ . The DM contribution is modelled as a 

'■J 

point source located at the central position of each dSph. 
In the analysis the free parameters include all the noma- 
lizations of the second LAT sources in the ROI, and the 
normalizations of the two diffuse backgrounds. 

Wc follow the method described in [3 to combine the 
results of different dSphs. Note that Jj factor is absorbed 



in when calculating . To get rid of the effect of 



Jj, we define 



fpi = Hj/Jj = 



87rm^ 



X Ci, 



(2) 



and write the joint likelihood function as 

j 

1 



In(lO) Jj^J^Tia 



- [logio ( -h ) -logio("'j)J 



(3) 



where D represents the binned 7-ray data and {p}j are 
the ROI-depcndent model parameters. log]^g(Jj) and aj 
arc the mean and standard deviations of the Gaussian 
distribution of logj^Q {Jj). The values of dSph J-factors 
and their uncertainties are given in Table I of 



http:/ /fcrmi. gsfc.nasa.gov/ssc/data 



^ http:/ /fermi.gsfc.nasa.gov/ssc/data/access/lat/Background- 
Models.html 
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FIG. 2: Comparison of the results using the likehhood map 
method developed in this work (solid) and that derived using 
Fermi Scientific Tool (long dashed). Here we assume 66 chan- 
nel of DM annihilation and the results from the ten dSphs are 
combined. The 2-year result from 0] is also shown. 
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We preform a grid scan with 10 energy bins logarith- 
mically spaced between 0.2 and 200 GeV, and 50 bins of 
i/'i logarithmically spaced from 10"^'^ to 10~^^ cm^^ s~^ 
GeV~'^. For each given (Ei, tpi) point, we vary all the 
nuisance parameters as well as Ji, known as the "profile 
likelihood" technique [lo| . and take the maximum likeli- 
hood of Eq. ([3]) as the final likelihood probability value. 
The combined likelihood map based on Ei and Efipi is 
shown in Fig. [TJ For each energy bin, we normalize the 
maximum likelihood to one. The solid line is a one-sided 
95% confidence level, corresponding to — 2Aln£ = 2.71. 

Combining the likelihood of each energy bin we can 
get the total likelihood for any input spectrum . The 
total likelihood can be calculated as 



(4) 



where tpi is the expected fiux of the spectrum in the ith 
energy bin, with w ^ . 

In Fig. [5] we give the comparison of our result with 
that derived using the Fermi tool (following directly the 
way in Q but using 4-yr data), for DM annihilation to 
bb channel. We can see that these two results agree well 
with each other. Also the two-year result given in Q 
is shown. The results show that the 4-yr constraint is 
even weaker than that of 2-yr constraint. Similar result 



was also reported recently 
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and it was possibly due 



to the update to pass 7 of the data and the statistical 
fiuctuations. 



FIG. 3: The 68% (inner) and 95% (outer) profile likehhood 
contour for the toy model. 



III. PARTICLE MODEL FIT 

In this section we apply the method described in Sec. 
II to specific DM models to derive the constraints. Three 
particle models, a toy model, the effective DM model, 
and the supersymmetry neutralino DM are adopted as 
examples. 



A. Toy model 

To generalize the upper limit on — (av) plane, 
we design a toy DM model which annihilates to 7-rays 
only via bb, t^t^, and W^W~ channels. We only con- 
sider these three channels because the differences between 
them are relatively large compared with other channels. 

With the 4-year likelihood map (Fig. [Ij, we con- 
duct a random scan of the four parameters, m^, {(jv)^i, 
{av)^+T.-, and (crw)vi/+iv- ■ The total cross section (av) 
equals to the sum of the three channels. The best-fit 
point obtained in the scan is: = 200 GeV, (av) = 
2.39 X 10~^^cm3s-\ BRf,b = 0.47, BR^+^- = 0.53, and 
BR^y+^y- = 10~^. Not surprisingly, the mixed chan- 
nels can give better fit due to more degrees of freedom. 
However, the current Fermi-LAT data can not effectively 
constrain the branching ratios, namely the shapes of 

In Fig. 131 we present 68% (inner) and 95% (outer) 
contours on the — (av) using the profile likelihood 
method. Since we project the model parameter space 
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FIG. 4: The 95% lower limits on A for each operator from the 
4-year Fermi-LAT dSph data. 



from 4 to 2 variables, the 68% and 95% contours are 



defined with /^L(D|7 



v)) dm^ d{av) equals to 0.68 



and 0.95, respectively. Here, R is the region within the 
contours. Compared with Fig. [2l we can clearly see that 
the 95% upper limit of this toy model is significantly 
higher than that of bb channel only for m^^ > 100 GcV 
and slightly higher for lower masses. This is because of 
the inclusion of r+r" channel in the toy model. Accord- 
ing to the results in for < 100 GeV the con- 
straints on {av)j£ and {av)T+^- are comparable, while 
for higher DM mass the constraint on {(jv)t+t- becomes 
much weaker than that on {av),,^. Including W~^W~ 
channel also makes the total constraint weaker. 



B. Effective DM models 

Before assuming particular DM models, we may adopt 
an effective interaction approach to describe t he i nte rac- 
tions of the DM particle with the SM particles 



12-17 . In 



such scenario, the DM particle exists in a hidden sector, 
which communicates with the SM sector via a heavy de- 
gree of freedom in the connector sector. At energy scale 
well below this heavy mediator we can introduce a new 
energy scales A to describe the effective couplings be- 
tween DM and SM fcrmions. In this study, we consider 
the interactions between fermionic DM field x ^-^d SM 
fcrmion / described by scalar (Ogcaiar), vector (Ovector) 



f ^ 



O.. 



scala 



E%^fe)(/7) , 



(5) 
(6) 
(7) 



where Ci^v,A,s is an effective coupling constant of or- 
der O that can be absorbed into A^ and / runs over all 
the lepton and quark states. Except ©scalar whose cou- 
pling also depends on to/, Ovoctor and Oaxiai both have 
universal coupling to all the fermions. 

We follow the calculation in [l^ to obtain (av) and 
for the effective DM models, and constrain each of 
the effective operators using the 4-year Fermi-LAT dSph 
data. Since we assume a universal coupling for all the 
fermions and the top quark channel only opens at > 
172 GcV, the upper limits of (av) of these three operators 
at low mass region uiy, < 172 GeV behave like bb channel 
in Fig. O Even if the ti channel opens, due to the large 
contribution to the 7-ray flux through the bb channel, the 
result is still not much different from the one xx directly 
to the bb. However, it is worthy to mention that the 
constraints on (av) of the three operators will depend on 
TO/ (sec Eqs. (Al), (A4) and (A7) in [3]). 

In Fig. 01 we show the 95% lower limits on {m^, A) 
plane for operators Oscaiar (black dash-dotted), Ovector 
(red sohd) and Oaxiai (blue dashed) and Ogcaiar- It is 
clear to see that both Oaxiai and Ogcaiar have a kink at 
the top quark mass because the tt channel opens. For 
operator Ovector there is only a very weak dependence on 
TO/, and the tt channel is suppressed. 

It is of interest to compare the limit on (m^ , A) plane 
driven by the Fermi-LAT gamma-ray low-latitude result 
[l6t and our dSphs lower limits. The limits from dSphs 
for all three operators arc slightly stronger. However, 
comparing our result with the one obtained from Galactic 
radio [17[, we found that the constraint for Ovector from 
dSphs can be stronger than the Galactic radio. 



C. The Minimal Supersymmetric Standard Model 

Undoubtedly, among various particle physics models, 
SUSY neutralino is the most popular DM candidate. 
However, many recent published results such Higgs 
boson candidate discovery at the LHC, flavour physics, 
(5((7 — 2)^ (the muon's anomalous magnetic moment), and 
the relic density of DM can set strong constraints on the 
SUSY parameter space. In addition, Fermi-LAT 7 ray 
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data can also test the SUSY DM, e.g. 
persymmetric standard model (MSSM) 18 

Sin 



the minimal su- 
19i and the 



Constrained MSSM 



In this subsection, we use 



4-year Fermi data on dSphs to test the MSSM scenario 
of DM. 

Given the large number of parameters in the MSSM, 
we have to take some reasonable and simplified assump- 
tions. We start by assuming that no CP violating phases 
are present, and that all mass matrices and triliuear cou- 
plings are diagonal in flavor space, in order not to vi- 
olate the quite stringent constraints on flavor changing 
neutral currents (FCNCs). Moreover, we take the first 
and second generations of scalars to be degenerate, again 
motivated by experimental constraints, such a.s K — K 
mixing. As the trilinears are proportional to the fermion 
masses, we can safely ignore the contributions from the 
first and second generations as they are dwarfed by the 
third generation. 

Unless there is a lighter gravitino or R-parity is not 
conserved, the lightest neutralino is the only MSSM par- 
ticle that can make a good DM candidate. The light- 
est neutralino is the lightest mass eigenstate of a mixed 
gauge-eigenstate, Bino, Wino, up-type Higssino, and 
down-type Higgsino, 



(8) 



The coefficients Zi (i = bino, wino, i/^) are deter- 
mined by diagonalizing the neutralino mass matrix. To 
describe the neutralino compositions, it is convenient to 
introduce a gaugino fraction, fg = .^^ino + ^wino- When 
fg is close to 1, gauginos will dominate the neutralino, 
on the other hand the neutralino will be higgsino-like if 

fs 

The input parameters and their prior ranges are shown 
in Table m After the Higgs candidate was discovered at 
the LHC 



22l . |23|, the Higgs resonance region, 2mp, 



ruh 



126 GeV, can drive a large 7-ray flux because 
of small and large {crv). Besides the Higgs reso- 
nance region, one can also have a large ipi at the Focus 
point region and A- funnel region with small m^. We 
let TOi = 0.5m2 to obtain a bino-like neutralino. How- 
ever, we still allow a Higgsino like neutralino, particu- 
larly at ~ 1 TcV, which is strongly disfavoured by 
5{g — 2)^ constraint. We keep the 1st and 2nd genera- 
tions of squarks as heavy as ^ 2.5 TeV and the first two 
generations of slepton masses mi = mi + 50 GcV at 
SUSY scale. 

To conduct the investigation, we use the nested sam- 
pling algorithm, implemented in MultiNest[24j which 
is incorporated into BayesFITS package, with 9000 live 
points, evidence tolerance factor 0.5, and sampling effi- 
ciency 0.8. To get rid of prior dependence and evenly 
explore the full parameter space, we make two separate 
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bino mass (GeV) 


mi = 0.5m2 


wmo mass (CjeVj 


10 < m2 < 4 X 10 


gluino mass (TeV) 


0.7 < m3 < 5 


top/r-quark trilinear (TeV) 


T ^ A A ^ '7 

~1 < At, Ar < 1 


r\ niiQvl^ "fT'ilmociT' ( ' 1 'cj \/ 1 
U LJLiaiJi LililliCcLi 1 _Lt! V J 


4,, — n 

jT-^ U.tJ 


pseudoscalar mass (TeV) 


0.2 < rriA < 4 


^ parameter (TeV) 


10^2 < M < 4 


3rd gen. squark mass (TeV) 


0.3 < < 4 


stau mass (TeV) 


0.1 < TTlf < 4 


lst/2nd gen. slepton mass (GeV) 


m? — mi + 50 

-1^1,2 


lst/2nd gen. squark mass (TeV) 


rrir, =2.5 

Ql,2 


ratio of Higgs doublet VEVs 


3 < tan ^ < 62 



TABLE I: The prior ranges of input parameters over which 
we perform the scan of the MSSM. We make two separate 
scans, one with flat priors and the other one with log priors 
for all mass parameters. In both scans, we adopt a flat prior 
for tan /3 and trilinear couplings. 

runs, one with fiat priors and the other with log priors 
for all the mass parameters. For the rest parameters we 
use only the flat priors. The observables included in the 
total likelihood >Ctotai used to constrain the parameter 
space are Higgs mass rrih ~ 126 GeV, WMAP relic den- 
sity n^h'^, Bs M^A* J b S7, Bu tv, and AMb^- 
For SM precision, we also include mw and sin^off- The 
theoretical and experimental errors can be found in Ref. 
25| and references therein. Regarding to the relic den- 
sity we assume that there is no other DM ingredients 
except neutralino. If the relic density is not dominant 
by neutralino, the DM fluxes must be rescaled by a fac- 
tor (ri;^/riwMAp)^ because of DM densities. We are only 
interested that £dSphs is improved by adding the rescal- 
ing factor but ipi is still large enough to be detectable 
by Fermi-LAT. To obtain the above solution, it is only 
possible that fi^ is slightly smaller than riwMAP- There- 
fore, in the level of scan, wc will only consider a Gaussian 
likelihood instead of an upper limit for neutralino relic 
density. In addition, we also check the exclusion bounds 
obtained from the Higgs searches at LEP and the Teva- 
tron, by implementing HiggsBounds-S . 8 . [l^. Note 
that we do not use 5{g—2)^ constraint because it strongly 
disfavoured for > 600 GeV where neutralino is most 
Higgsino-like. 

The total 7-ray spectrum from neutralino annihilation 
is given by 



f 



X 



dE 



(9) 



where X runs over the contributions of secondary 
photons, final state radiation, virtual internal 
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FIG. 5: The impact of 4-year Fermi-LAT dSph result on 
(m^, (av)) plane of neutralino DM. The grey dots satisfy 
n/i^ < 0.224 and — 21n£totai < 20, and the blue squares 
are selected from grey cycle points with rrih > 120 GeV and 
B < 10. 
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FIG. 6: The scatter plot on the (m^, a^^) plane in the MSSM. 
The color scheme is same as Fig. [5] The solid (red) line 
is XENONIOO (20f2) 95% upper limit and the dashed (ma- 
genta) line is XENONfT projected sensitivity. 



bremsstrahlung, and monochromatic 7-ray line, Bf 
is the branching ratio into the annihilation channel /. 
The contribution from secondary photons is mainly 
produced through the decay of 7r° and K°. For com- 
pleteness, we also include the monochromatic 7-ray 
line contributions because XX ~^ 11 1 1^ can slightly 
contribute to our likelihood for niy^ < 200 GeV. The 
spectrum is calculated using the DarkSUSY-5 .0.6 



package [27 1. 

It is possible that the annihilation cross section of 
DM in the Milky Way today can be boosted by non- 
perturbative Sommerfeld corrections compared with that 
in the early Universe, especially for > TeV in the 
MSSM case (e.g., HB-Hg). In the MSSM, the light- 
est neutralino with mass ~TeV arc mostly Wino- or 
Iliggsino-like, with Sommerfeld enhancement factor of 
order 0(1) 28-3^. On the other hand, can dif- 
fer by a factor of several between different simulation 
codes PythiaS [31|, PythiaS [32| and Herwig [33|. To 
investigate possible source of boost factor which could 
come from either the Sommerfeld enhancement or un- 
certainties of we can employ a phenomenologi- 
cal boost factor to boost our SUSY prediction flux 
,0SUSY order to obtain the maximum dSphs likelihood. 

We present the scan result on (m^, (cru)) plane in 



Fig. [5] Totally ^ 200, 000 sample points are obtained. 
The grey dots correspond to the sample surviving from 
-21n/:totai < 20 and ^h^ < 0.224 cut. From Ref. ji^, 
we can see Sommerfeld corrections can also reduce the 
relic density. Since we do not include the Sommerfeld cor- 
rections in our relic density computation, we here apply 
a broader range flh^ < 2f2wMAP^^ cut in order to con- 
servatively include the possible points. The blue squares 
arc selected from the grey dots with criteria > 120 
and B < 10. We find that the blue squares for m^^ < 300 
GeV are mostly bino-like neutralino but for to^ > 1200 
GeV they are mostly Higgsino-like neutralino. For those 
Higgsino-like neutralino the boost factor is found to be 
S ~ 4, which could be due to either the Sommerfeld effect 
or the uncertainties of For the bino-like neutralino 
with lower masses, few points are found with S < 4. 
Since the Sommerfeld effect does not apply on such small 
DM mass range and the uncertainty of is also only 
a factor of several, we conclude that the 4-year Fermi- 
LAT data is less sensitive to bino-like than Higgsino-like 
neutralino. 

The constraint on the spin-independent cross section 
between DM particle and nucleon a^^ is shown in Fig. 
El Also shown are the XENONIOO (2012) result of the 
direct detection experiment [3J] and the expected sen- 
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sitivity of XENONIT. We notice that XENONIOO can 
exclude most of neutralinos which have gaugino fraction, 
0.6 < fg < 0.9. However, for pure gaugino ~ 1 or Hig- 
gsino fg ~ 0, Fermi-LAT data still can test them prior 
to XENONIT. 



IV. SUMMARY 

In this work we develop a new method to ensure fast 
computation of the likelihood of any spectral component 
from Fermi-LAT observations of dSphs. The Fermi-LAT 
data of each ROI are binned into several energy bins and 
the likelihood of a point source contribution located at 
the ROI center in each energy bin is calculated with the 
public Fermi Scientific Tool. A likelihood map which de- 
pends only on the Fermi-LAT data can be derived with- 
out the speculated energy spectrum of the new compo- 
nent (Fig. [I}. Combining the likelihood of each energy 
bin we can get the total likelihood of any given energy 
spectrum. This method is tested to be in good agree- 
ment with the results directly derived with Fermi Scien- 
tific Tool. It is thought to be of great convenient for the 
search for DM signals from the Fermi-LAT data, espe- 
cially for the model-dependent search. 

We apply such a method on several DM models, in- 
cluding a toy model with mixture of several channels, 
the effective operator scenarios and the MSSM models. 



The Fermi-LAT data can give interesting constraints on 
the DM models. First, in our toy model fit, we found 
that the current Fermi-LAT data cannot effectively con- 
strain the shapes of Second, the Fermi-LAT data 
can improve the lower limits of three effective operators 
on (m^,A) plane. Finally, for the MSSM, we find the 
Fermi-LAT data are more sensitive for the Higgsino-like 
neutralinos than the bino-like neutralinos. Furthermore, 
Fermi-LAT data can even test the MSSM models with 
pure gaugino or Higgsino compositions more effectively 
than the current direct detection experiments. 

The data of the likelihood map and a FORTRAN code 
35j to calculate the final likelihood of any input spectral 
function arc provided for download. 
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